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1. Introduction 

The famous Black-Scholes model with its relatively stringent assumptions does 
not capture many phenomena of modern financial markets. A prime example is 
the stochastic nature of the financial asset's volatility, called volatility smile (see 
for example Hull and White [6]). In recent years many stochastic volatility models 
have been introduced and developed, but as was to be expected, making the volatil- 
ity stochastic has complicated the models considerably (see for example Rebonato 
|16j). It is not our aim to review the broad range of stochastic volatility models but 
rather to focus on and develop the idea of modeling stochastic volatility in the sim- 
plest possible but effective way. SABR is an excellent example of a model complex 
in nature but simple in form. This well known and celebrated model, introduced 
in 2002 by Hagan et al. [5j, has since been effectively used and investigated by 
market practitioners. Soon after its introduction it turned out to be more effec- 
tive than Black-Scholes and local volatility models. The key idea in SABR is to 
make stochastic volatility a simple stochastic process and then shift the difficulty 
of finding the financial asset's distribution to the level of finding the distribution of 
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the diffusion describing the asset price. For SABR, determining closed formulae for 
the asset distribution remains, in general, an unsolved task (as far as the authors 
know) . 

The task of determining closed formulae for the probability distribution for 
SABR with parameter beta equal to one, called a log-normal stochastic volatil- 
ity model, has recently been investigated by Maghsoodi [H], [12] ■ In this case it is 
not particularly difficult to write out the solution of the model, i.e. the stochastic 
process representing the asset price, as the exponential of a linear combination of a 
pair of correlated Brownian motion functionals. Maghsoodi used the techniques of 
changing time and measure to find the joint density function of these functionals. 
The same techniques had been used earlier by Yor in the problem of valuation of 
Asian options (see (25). However, Maghsoodi has not mentioned that the asset 
price loses the martingale property in a log-normal stochastic volatility model in 
the case of positive correlation between the asset and its volatility. 

In our work we propose a kind of reversing the idea of the SABR model and 
continue the line of research of Hull and White [6] followed also by Romano and 
Touzi [18] as well as by Leblanc [lO]. We shift the complicated nature of the model 
to the level of the process representing volatility, keeping the diffusion of the asset 
relatively simple. So, we assume that the asset price process X satisfies dXt = 
YtXtdWt with Y given by dYt = ^{t, Yt)dt + a{t, Yt)dZt, where the processes W, Z 
are correlated Brownian motions. We call this model a linear stochastic volatility 
model. Under some natural assumptions we are able to prove that the distribution 
of the asset has a density function and derive the probabilistic representation of that 
function (Theorem 12. 2|) . This representation depends on some functionals of the 
process representing volatility, so the problem of determining the asset distribution 
reduces to finding the distribution of a 2-dimensional functional of the volatility. 

In Section 3, we point out two nontrivial examples of such models in which we can 
benefit from probabilistic representations of the asset density function. The first 
example is a log-normal stochastic volatility model which is a SABR model with 
beta equal to one. Using the result of Matsumuoto and Yor [14] who derived the 
density function for the vector of Brownian motion with drift and its exponential 
functional, we find closed formulae for the density function in a stochastic log- 
normal volatility model. The second example we present is the one with volatility 
being a 3-dimensional Bessel process (briefiy, a BES(3)) starting from 1. Following 
the ideas of Pal and Protter [15] as well as Donati-Martin and Yor [3], we use 
change of measure techniques and obtain the Laplace transforms for the volatility 
functionals. Next, we derive closed formulae for the density function in this model, 
using Theorem 12.21 

In Section 4 we derive probabilistic representations for European call and put 
option prices in linear stochastic volatility models. The probabilistic representation 
for vanilla option prices is independent of the distribution of the asset price itself. 
This allows us to obtain formulae for the arbitrage prices of vanilla options in a 
log-normal stochastic volatility model. 

Similar probabilistic representations for European call and put option arbitrage 
prices in a linear stochastic volatility model have also been given by Romano and 
Touzi [18], but in a slightly different context. They considered a slightly different 
model and established a set of assumptions under which they obtained probabilistic 
representation results while proving the convexity of European call and put options 



REPRESENTATIONS IN A LINEAR STOCHASTIC VOLATILITY MODEL 



3 



in their setting (also linear in our sense). In particular, they assumed that the 
coefBcients ^ and a in the definition of Y are bounded. In our work we relax this 
assumption (see Theorem 14. ip . In our examples the drift coefficient is not bounded, 
but the probabilistic representation for option prices holds. It is easy to check that 
in the case of the volatility process Y being BES(3) with Yq ~ 1 the condition 
p < implies that X is a martingale, and the probabilistic representation works 
well in that case. The same is true in the case of a log-normal stochastic volatility 
model. Thus for linear stochastic volatility models our result is more general than 
the one in [18]. It should be mentioned that Leblanc [10] gives the arbitrage price 
of call option in a linear stochastic volatility model, with some concrete examples 
of volatility, in terms of Laplace and Fourier transforms. 

Closed formulae for the density function and vanilla option prices in a stochastic 
log-normal volatility model are interesting and important for appHcations since 
such models are popular, especially among forex exchange options traders (see [5]). 
Similar results for log-normal stochastic volatility models were also presented in [11] 
and [12] . It should be mentioned that during the writing of this text we discovered 
that Forde [4] also found the link between the result of Matsumuoto and Yor [14] 
and SABR for beta equal to one, but he did it only for uncorrelated price and 
volatility. Thus our result for a log-normal stochastic volatility model is deeper 
than that in [4] as it requires only the asset price be a martingale. 

In the Appendix we give a conditional version of the Donati-Martin and Yor 
formula [^ and establish the density function of the vector (i?t,/p Bf^du), where B 
is a standard Brownian motion. 

2. PROBABILISTIC REPRESENTATION OF THE DENSITY FUNCTION OF THE ASSET 
PRICE IN A LINEAR STOCHASTIC VOLATILITY MODEL 

2.1. A linear stochastic volatility model. We consider a market defined on 
a complete probability space (ri,JF,P) with filtration F = {^t)t£[o,T\i T < oo, 
satisfying the usual conditions and T = Tt- Without loss of generality we assume 
the savings account to be constant and identically equal to one. Moreover, we 
assume that the price Xt of the underlying asset at time t has a stochastic volatility 
Yt, and the dynamics of the vector (X^Y) is given by 

(1) dXt = YtXtdWu 

(2) dYt = fi{t, Yt)dt + a{t, Yt)dZt, 

where Xq, Yq are positive constants, the processes W, Z are correlated Brownian 
motions, d{W, Z)^ = pdt with p € (-1, 1), and p : M+ x E+ ^ M, cr : M+ x M+ ^ M 
are continuous functions such that there exists a unique strong solution of ([2|, 
which is positive and Y^du < oo P-a.s. 

Under these assumptions the process X has the form 

(3) Xt=Xoe/o>'"''i^-/o'>'.^rf«/2^ 

and this is a unique strong solution of SDE ([TJ on [0, T]. The existence and unique- 
ness follow directly from the assumptions on Yt and the well known properties 
of stochastic exponent (see, e.g., Revuz and Yor [E]). The process X is a local 
martingale, so there is no arbitrage on the market so defined. 
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We call this model a linear stochastic volatility model, because the SDE ([T]) 
governing the asset price is linear with respect to the asset itself with coefficient 
being the stochastic volatility Y. 

Remark 2.1. a) It is worth mentioning that the constant p in the model can be 
replaced by a measurable, deterministic function p : [0, T] (—1, 1) and the results 
of this work remain true with minor modifications. 

b) Our standing assumption is \p\ < 1. However, our methods allow finding the 
distribution of Xt in the case p = ±1 (see Remark [53]). 



2.2. Existence of the density function and its probabilistic representation. 

We start with the main theorem of the paper on existence of the density function 
of the underlying asset in a linear stochastic volatility model, and its probabilistic 
representation. This representation allows us to find a closed formula for the density 
function (see examples in the next section), which is important for applications (see, 
e.g., Carmona and Durrleman [2]). 

Theorem 2.2. Fix t <E [0,T]. // 

(4) E / Y-'^du < oo, 



then the random variable Xt has a continuous density function gxt ■ Moreover, the 
density function has the probabilistic representation 

1 .^'(^^~^^^*^ 



(5) 9xAr)=^ ... ^ , ... 

[raz(t) V o-^l*) 

where <f> is the cumulative distribution function of a standard Gaussian random 
variable N{0, 1), and 



(6) pzit) = p I YudZu - \ / Y^du, 



(7) al{t) = il-p') fY^du. 

Jo 

Proof. Notice that we can represent W in the form 



(8) Wt = pZt + - P^Bt, 

where {B, Z) is the standard two-dimensional Wiener process. 

Step 1. Proof of the existence of density. We prove that the characteristic function 
(j)Ut of Ut := \a.Xt is integrable. The Ito lemma applied to |(T]) together with ([2]) 
and |[8]) implies that 

\nXt^\nXf, + ez{t) + eB{t), 

where 

ez{t) ■.= p f YudZu - Ip^ f Y^du, 
Jo ^ Jo 

dsit) := v/r~^ /* Y^dB., - i(l - p2) r 
Jo ^ Jo 

Let = a{Zu ■.u<t). For s e R, 



du. 



\(l)uAs)\ = |E(e*''''^(*)E[e"^''«(*)|J^f])| < E|E[e*'*''«(*)|J^j 
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Since SDE ^ has the unique strong solution, there exists an appropriately mea- 
surable function ^{ , ) such that Y = 'i>{Yo,Z) (see e.g. |7]). Together with the 
fact that the processes B and Z are independent Brownian motions, this implies 
that the random variable OB{t), for a fixed trajectory of Zu, u < t, has Gaussian 
distribution with mean 



and variance 

Consequently, P-a.s. 
and hence 



4" = (!-/.») ['y, 
Jo 



■du 



du. 



E 



Therefore, by the Fubini theorem, we obtain 



Ac/,(s)|d.s< 



= E. 



27r 



(l-p2) 



< 



27r 



(l-p2 



In the last inequality we have used the Jensen inequality for the concave function 
/(x) = yjx. Now, using the integral Jensen inequality for the convex function 
/(x) = for a; > 0, we have 

ft 



1 



(9) 

Hence and by assumption (0]), 



< ^ / Y-'^du. 



27r 



E 



< 



27r 1 



(l-p2)t 



-WE / Yu du < oo. 



In such a way we have proved that the characteristic function (j)Ut is integrable. 
By the Inverse Fourier Transform theorem there exists a bounded and continuous 
density function of Ut- Since Xt = e^^ , we conclude that there exists a continuous 
density function of Xt . 

Step 2. Proof of ^ — probabilistic representation of Xt. For fixed r > 0, using the 
same arguments as in Step 1 we have, by ([3]) and |[8]), 



fiXt < r) = El r / , rt . \ \ 



EE 



(10) 



^niiz{t) + Gz{i)g<\^ — \Tt 



m[ - 



Xo^ 
h{rl 



EP g < 



<Jz{t) 
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where g is a standard Gaussian random variable independent of Tf , fJ-z{t) and 
<T|(t) are given by ^ and (O, respectively. The function h, defined for r > by 
(fTOl) . is continuous. Moreover, 



dr 



< 



(7z{t) 
1 



raz{t) 



cTz{t) 



2Trraz{t) 



< 



27rr 



27rr 



where 



27rr 



{l~p')t' Jo 

In the last inequality we use ((9]). From assumption |[4]) it follows that the random 
variable /i is integrable. Consequently, 



(11) 



Kir) - |-P(Xi < r) 



1 ^^/l^il-^^W 



razit) 



Now we prove that K is a density. It suffices to show that h[{r)dr = 1. Using 
the Lebesgue dominated convergence theorem we get 



(12) 



limE$ - 



In^-Mz(t) 



Hence, by (fTTjl and the Fubini theorem, 

^n^ 



E 



r-»o V az{t) 



0. 



.^0 



1 



l^zit)^ 



E 



dr 



10 u, X az{t) 
so gxt{f) ~ ^t(^) is the density function of Xf. 



E 



dr 



^/ In^-Atz(t) - 







) 


- 



□ 



Remark 2.3. From the last theorem it is clear that finding the distribution of Xt, 
for fixed i, reduces to deriving the distribution of the vector (/J YudZ^, Y^du), 
provided 10]) is satisfied. 

Remark 2.4. Note that in the case p = ±1 we cannot write the formula ^ for 
the density of Xt. However, we have W = ±Z in this case and 



(13) 



Xt 



SO the problem of finding the distribution of Xt, for fixed t, reduces, as previously, 
to deriving the distribution of the vector (/^ YudZu, Y^du). In the case of a 
lognormal stochastic volatility model we can use the results of Matsumoto and 
Yor [14] to obtain the distribution of [J^ Y^dZu, Jg Y^du), as we can express its 
components in terms of At and Vt just as in the proof of Theorem 13. II and use (fT9]) . 
In the case of volatility being BES(3), Theorem 13.31 gives us the Laplace transform 
of the vector {J* K„c?Z„, Y^du), so after inverting it we obtain its distribution. 
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3. Application of the probabilistic representation 

Using the probabilistic representation of the density, we can find a closed form 
of the density function in some cases including trivially the Black-Scholes model, 
and nontrivially a log-normal stochastic volatility model and a linear stochastic 
volatility model with volatility being a 3-dimensional Bessel process. 



3.1. Closed form of the density function in Black-Scholes and log-normal 
stochastic volatility models. The Black-Scholes model and log-normal stochas- 
tic volatility models are special cases of linear stochastic volatility models. A log- 
normal model was considered by Hull and White in the case of uncorrelated noises 
[6j, and it is a SABR model with /3 = 1, introduced in 2002 by Hagan et al. [5], in 
the case of correlated noises. 

The case of the Black-Scholes model is trivial. The process Yt = a > represents 
the volatility of the asset, p = 0, so it is a linear stochastic volatility model and 
Theorem 12.21 gives the well known density function of a random variable with log- 
normal distribution. 

A log-normal stochastic volatility model is much less trivial. In this case the 
functions appearing in the SDE for volatility are ^{y) =0 and a{y) = cry for y > 0, 
where cr is a positive constant. Thus the process y is a geometric Brownian motion 
and 



(14) 



Yt = Foe 



aZt~cr^t/2 



Moreover, it is easy to check that for any t <T, 



E / Yfdu 



^0 ^2 



1] < oo, 



so a log-normal stochastic volatility model belongs to the class of linear stochastic 
volatility models. 

Our main goal in this subsection is to find, for a log-normal stochastic volatility 
model, a closed form of the density function of the random variable Xt for fixed 
nonnegative t (see [12] for another result in this direction). We determine the true 
distribution of the price process, so we find a simple way to price derivatives in that 
model. 



Theorem 3.1. In a log-normal stochastic vo 
the price Xt of the underlying asset has the form 



model the density function of 



9Xt (r) 



where 

(15) 

(16) 



OO poo 



1 ^^,n^Y-.-fi^^y) + Yoy^ 



— s— 



f{x,y)^!-Yo[e^-l]-£^Y,fy, 



Gt{x,y) = cxp 



X t 1 + e^^ 
2^8 2y~ 



Gta^-{x,y)dydx, 



y J y 
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and the function 9 is defined, using hyperbolic functions, by the formula 



(17) e{r,t) 



\/2ttH 



-4^/2t-rcosh(^) 



sinh(^) sin 



Proof. Set Yt ■— Yt/a'^- It is clear, from fT4|) . that 

where Zt = aZtj^i is a Brownian motion. We can express and o\{t), defined 

by ^ and in terms of %: 



mit) - -Yo]-^ I Y^du, alit) = 



1-p 



2 rto- 



Y.du. 



Let 



14 



Then Yt = Yoe^' and /g Y^du = Fg^^t. Since 

1 , 



E / Y-'du ^ 



iYoa^ 



^3a-t _ 1) < 



we can use Theorem 12.21 and write the density function <?Xj, 2 in terms of Vt and 



At: 



(18) gx„.{r)=E 



rYo^/A^ 



-^0 



f{Vt,At) + Y,^At'-§^ 



YoJAt'-f 



where / is given by l|15p . Now, we use the result of Matsumoto and Yor |14| which 
gives the density function of the vector {Vt,At): they proved that for t > 0, y > 

and a: e M, 



(19) 
where 



P{Vt e dx,At e dy) = Gt{x,y)dxdy, 



Gt{x,y) = exp 



X t 1 + e^^ 
2^8" 2y 



— ,t -, 
V J V 



sh(0 



V27r3i 

Hence (fT8l) can be written in the form 



sinh(^) sin — I d^. 



00 nOO 



00 JO 



1 _^,n^i-o-f^^^y) + ny'-^ 



Gt{x,y)dydx, 



with /, G given by (fT5|) and (fT6|) . Replacing t by icr^ in the above formula finishes 
the proof. □ 

Remark 3.2. Although the formula for the density function of the price in the 
log-normal stochastic volatility model is complicated, this result describes the true, 
not approximate, probabilistic law for Xt. If X is a martingale, so describes the 
arbitrage price of the asset, having the density function we are able to use the 
risk-neutral valuation formula to price attainable European contingent claims. For 



REPRESENTATIONS IN A LINEAR STOCHASTIC VOLATILITY MODEL 



9 



example, evaluating the arbitrage price of power option (see, e.g., [19]) reduces, by 
Theorem 13. H to computing the integral 



with /, G given by l|T5|) and l|T6|) . We stress that in this way we reduce the valuation 
problem to numerical integration of the derived density function, as is usual in the 
literature (see e.g. [2]). Thus we avoid using asymptotic expansions (as in [5]); 
however, some difficulties arise during the numerical integration (see e.g. [l]). They 
are caused by the oscillating nature of the so called Hartman- Watson distribution 
density function which is a part of the density function derived by Matsumoto and 
Yor [11]. 

3.2. Closed form of the density function in the case of stochastic volatility 
being a 3-dimensional Bessel process. We now consider another example of a 
linear stochastic volatility model. We assume that the process Y, representing the 
volatility of the financial asset, is a 3-dimensional Bessel process with Yq = 1. It is 
well known (see, e.g., Revuz and Yor [17]) that F is a strong solution of the SDE 

(20) dYt = dZt + ^dt, 

where Z is a one-dimensional standard Wiener process. Since Yq = 1 and for any 

t < T, 



(21) E 



[ Y^du =/"(! + ?,u)du = < + < oo, 
Jo Jo 2 



the SDE's (dl) and i(20)) define a Hnear stochastic volatility model (with fi{y) ~ l/y 
and a{fj,) = 1 for y > 0). Moreover, it is easy to check, using the ltd lemma, that 
condition ^ is satisfied: 

ft 1 

(22) E / — du^ElnY^^ < KY^^ ^l + M<oo. 

Jo 

Therefore, as we noticed in Remark [231 the problem of finding the density function 
of the asset reduces to deriving the distribution of the vector (/p YudZu, Jq Y^du). 
We do this in the next theorem finding the Laplace transform of this vector. 

Theorem 3.3. Let Y be a BES(3) process with Yq ~ 1. The Laplace transform of 

the vector {J* YudZ^, /g Y^du) is given by the formula 

(23) 

Eexp(^-A YudZu-- Y^duj = p{t,y)yexpi^-X—^-^jmAy)dy, 



where A > 0, 6 > 0, 

p{t,y) = 



'2Tlt 



'^"Pl -^^j-""PV 2t 
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and the function H^^t is defined for y > by the formula 
(24) Hb.t{y) = 



bt 



exp ( - l{(y^[6tcoth(M) - 1] + 2btiy + D^^^^g^ 



sinh(M) 



Proof Note that the Laplace transform is well defined, since /* Y^du > and 



(25) 



ft I 

Y^dZu = -[Y^-l-2,t]>- 



1 + 3T 



Take a probability space (ri,jr, Q) and a process Y which is, under Q, a Brownian 
motion starting from 1. Define tq = 'mi{t > : Yj = 0} and let P be a new 
probability measure given by 



Yt 



t/\To- 



Then, by the Girsanov theorem, Yt is BES(3) with Yq = 1 under P (see Pal and 
Protter [13). Fix A > and 6 > 0. Then, by ^ and the fact that FtAro = for 

t > To, 



h2 rt 



E(^cxp I - A / YudZ, 
(26) = E, 



2 7o ^""^^ 



YtAro exp - A y„ [dYu - Y-^'rJu] - ^ £ Y^du^ 

le measure P, a Broi 
inf{t > : Bt = 0}. Then §^ impHes 



Let B be, under the measure P, a Brownian motion starting from 1, and set Tq 



u2 rt 



E exp -A / YudZu- 



E 



Y^du 



Bf^^B exp ( — A / BudBu + Xt — 



h'2 rt 



Btdu 



^ Jo 



If T(f < t, then the expression under expectation vanishes as B^b — 0. Hence we 
can replace i by t A r^f under expectation, so 
rt i2 rt 



E 



B^^^B exp - A / B„dB^ + At - — / B^rfu 



E 



^tAr„« exp - A 



(27) 



E 



^tAr,? exp 



A 



t^ri-l^Bl^. 



E 



2 /•tAr„'^ 



exp 1 - y 



B, 



From Theorem 15.11 (see Appendix) for a; ~ 1 , it follows that 

62 rt 



(28) 



E 



exp 



B?,du 



Bt = y 
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where Hb^t is given by l|24|). Inserting l(28|) into ((27l) we obtain 



E( exp( - Ay y„dZ„ -y^ 



Y,fdu 



t A Tr 







^ -j^ _|_ ^2 

l{ro«>t}l{St>0}5tCXp ( - A ^ ]Hb^t{Bt) 



J l,y)yexp ^ - A 



i - 1 + 2/' 



Hb,tiy)dy, 



where 



exp 



(y-1) 
2t 



exp 



2i 



is the transition function of Brownian motion starting at 1 and absorbed at zero 
(see Karatzas and Shreve [U p. 97]). This concludes the proof. □ 

The last theorem provides the Laplace transform of the vector (/J YudZu, /q Y^du), 
so we are able to invert it and obtain the density function in a linear stochastic 
volatility model in the case of the stochastic volatility being a 3-dimensional Bessel 
process. However, as shown below, we can solve this problem in another way. The 
difference between the two methods of finding the density function is that in the 
second one we do not need to invert the Laplace transform of a vector but of a 
function of one variable. 

We start from a lemma computing EF{Yf^,J^^ Y^du) for a measurable function 
F and for Y being EES (3). 

Lemma 3.4. LetYt be a BES(3) process starting from 1, and g be the density func- 
tion of the vector {Bt, B^du) , where Bt is a standard one- dimensional Brownian 
motion. Then for a measurable function F such that KF {Y^ , Y^ du) < oo, we 
have 



OO rOO />OC nOC nOO pOO 



OO J J —oo J 



Q*ixi,yi,X2,y2, X3, ys)dxidyidx2dy2dxzdy3 



whe 



''{xi,yi,X2,y2,X3,y3) = F{1 + x* ,t + y*)g{xi,yi)g{x2,y2)g{x3,y3), 



and X* = xl + xl + x\, y* ^yi+y2 + ys- 

Proof. If Yt is a EES (3) process starting from 1, then 



Y,^^l + Zf, + Zl, + Zl„ 
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where Zi^f , Z2,t, ^3,t are independent standard Brownian motions (see [HI)- The 
conclusion of the lemma now follows from the chain rule for conditional expectation: 



m.F{\ + zl, + zl, + zl,, t + j\zl, + zl^ + zljdr^ i^f-^^ 

/ g{xi,vi)^F[l + xl + Zl^ + Zl^, t + yi+ [Zl^ + ZlJdu)dyidx 

-ooJo V Jo J 



□ 



Using this lemma we find the density function of the asset price Xt. 



Theorem 3.5. Let Yt be a BES(3) process with Yq = 1, and g he the density func- 
tion of the vector {Bt, B'^du), where Bt is a standard one- dimensional Brownian 
motion. Then the density function of Xt, for fixed t >0, is given by 

/OO pOO pOC nOC pOC pOO 
/ / / / / Qir,Xo,xi,yi,X2,y2,X3,y3)dxidyidx2dy2dx3dy3, 
-ooJQ J-ooJO J-ocJO 

with 

Q{r,Xo,xi,yi,X2,y2,X3,y3) 



- 3i] - ^[t + y*]\ g{xi,yi)g{x2,y2)gix3,y3) 



^il-p^)[t + y*] J r^{l-p^-)it + y*) 

where x* = x\ + x\ + x\, y* = yi + y2 + 2/3 and $ is the cumulative distribution 
function of a standard Gaussian random variable N{0, 1). 



Proof. Since (|22|) is satisfied, we conclude from Theorem 12.21 that 



1 



-r^{l-p^)jlY^du V ^{l-p^)jlYidu 
and the assertion follows immediately from Lemma [ 



□ 

Remark 3.6. In order to complete the example of BES(3) volatility (starting 
from 1) we comment on the density function g of the vector [Bt, B^^du), where 
Bt is a standard Brownian motion. It is proved in the Appendix (see Corollary 15. 2p 
that for c> 0, i > 0, X e R, 



(29) 
where 

(30) 



f 

Jo 



l{x,y)dy = H*(t,c, a;), 



H*{t, c, x) 



1 / V2c / a;2 



^ cxp \/2ccoth(i\/2c) 

2tt V sinh(^^/2^) V 2 ^ ^ 



So if we denote by £ ^ the inverse Laplace transform, then 
(31) g = C-'H* 
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and we obtain the density function of the vector (St, /q B'^du) 



4. Closed form of the arbitrage price of a vanilla option in a linear 

STOCHASTIC volatility MODEL 

In this section we derive a probabiHstic representation of a vanilla option price 
in a linear stochastic volatility model. We are interested in computation of the 
arbitrage prices, so the process X describing the discounted price of the asset should 
be a martingale. Next, as examples, we show how to deduce from Theorem 14.11 
closed formulae for option prices for the models of Section 3. In our examples we 
give conditions guaranteeing that X is a martingale. Then, just as in Section 2, we 
show how the valuation of vanilla options in that model can be reduced to finding 
the distribution of the vector (/^ YudZu, Jq Y^du). 

4.1. Probabilistic representation of the arbitrage price of a vanilla option 
in a linear stochastic volatility model. Under the assumptions of Theorem 12. 21 
and for the process X being a martingale, we provide probabilistic representations 
for the arbitrage prices of European call and put options. These formulae generalize 
the famous Black-Scholes formulae as well as the result of Hull and White for a 
stochastic volatility model with uncorrelated noises [^. 

Theorem 4.1. Assume that X is a martingale. In a linear stochastic volatility 
model, under assumption ^ , the time zero arbitrage prices of European call and 
put options with strike K > and maturity t have the following probabilistic repre- 
sentation: 



(32) E[Xt-K]' 



Inf 



In^ 



-Mz(i) 



(33) nK~Xt] + 

fhi^ - lJizit)\ \ ,,.^2,,.,^ f\n^- iiz{t)-al{t) 

V ^z{t) ; M V ^z{t) 

where nzii) and <7^{t) are given by ^ and 

Proof. Let K > 0. We compute E[K — Xt]^ . Integrating by parts, using the 
probabilistic representation of gxt, i.e. ([5]), we have 



E[K-Xt]+^ J {K -r)+gxAr)dr = J (/f - r)+— E$ j^- 



az{t) 



dr 



/In^-fizit) 



K 



K 



10 V az{t) 

where in the last expression we have used (fT2l) . Next, by the Fubini theorem 



E$( — ^" Z^^"" ' \dr = E 



/o 



^z{t) 



fir ~| dr 

\ r>Xoe"'zW+t'zW ^ 



ds = h, 
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where 



Let us define 



K* 



K 



-tJ.z(t) 



ds. 



In 



+ /iz(t)+a|(i) 



In 



K 



Then, changing variable s to v ^ s(Jz{t) in /s, we obtain l(33|) : 

1 



/3 =E(^Xoe^^(*) 



— OO V 2naz(t) 
\nK* 



g4(t)/2$ 



A'* 



Hence, as X is by assumption a martingale. 



InX* - cr|(t) 



(34) 



C[Xt - K]+ =Xa- K - E[/f - Xt 



= Xo(l-E[, 



,Mz(t)+fT|(t)/2 



z(*)/2$(_d*)]) - K(l - E$(-d^)). 



Let us observe that 



Moreover, 

(35) Ee/o^"rfW'"-/o>'.'rf«/2 



iE 



E 



gP/„*y„rfz„-pVn*y„^du/2^ 



Since X given by ([T]) is a martingale, we have 1 = Ee-'^t) ^^udw,, /(, y^d«/2 ^^j^j ^j^g 
process (e^^(*)+°'2(*)/^)t is a martingale. Hence by (|34| we finally conclude that 



E[Xt - Ky- = XoE|^e''^^'^+''2^'^/' (^1 - $(-dD j J - - ^(-^2) 

= XoE e^^(*)+'"^(*)/2$(d*) -ifE$(d;), 

which ends the proof. □ 

4.2. Examples. In this subsection we consider the previously discussed models. 

4.2.1. Black-Scholes and stochastic log-normal volatility models. In these two cases, 
closed formulae for the arbitrage price of European call and put options with strike 
K > can be derived. We emphasize that these results are not a direct consequence 
of deriving the density function for the model. Rather, they are consequences of 
the probabilistic representation (see Theorem 14. ip of the arbitrage price of vanilla 
option in a linear stochastic volatility model. 

In the case of the Black-Scholes model, nz{t) ~ — and cr|(t) = cr^t, so ([32|) 
and l(33|) immediately give the famous Black-Scholes formulae. 
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As before, the case of a log-normal stochastic volatility model is less trivial. We 
give formulae for the arbitrage prices of vanilla options in such models (different 
formulae were obtained in [12] in another way). 

Remark 4.2. .Jourdain [8] proved that the condition p S (—1,0] is equivalent 
to X being a martingale. So, in further considerations, whenever we need X to 
be martingale, we consider only nonpositive p, and in this case P is a martingale 
measure. 



Theorem 4.3. In a log-normal stochastic volatility model with p < 0, under as- 
sumption ^ , the time zero arbitrage prices of European call and put options with 
strike K > and maturity t are given by 



(36) E[Xt-K]- 



(37) E[K^Xt]- 



OO nOO 



OO Jo 



Xoe-^(^'«)$(di(x, y)) - K<^{d2{x, y)) 



OO poo 



oo Jo 



K<i>{-d2{x, y)) - Xoe^(^'^)$(-di(x, y)) 



Gta^{x,y)dydx, 



Gt^2{x,y)dydx, 



where f, G are given by lfT5|) and (fT6|) . and 



lnj^ + fix,y) ^ Yo 1 - 
di [X, y) = ; ^ „ + -TT \/ —y, 



d2{x,y) = di{x,y) - / y. 



Proof. The assumption p < and Remark 14.21 implv that X is a martingale. Ar- 
guing as in the proof of Theorem 13.11 and using the same notation we have, by 
Theorem 14.11 



(38) nX^-K]+=E[Xoef^''^'^*^^{d^{Vt,At))~K<^{d2{VuAt))\, 

(39) E[A' - X^]+ = E[ - K<^>{-d2{Vt,At)) - Xoe^(^-^')<I>(-di(F*, ^t))] , 

and hence 
(40) 



OO pOO 



- K] 

(41) 

/OO /"OO 
/ 
-oo J a 

To conclude the proof we replace t by tcr^ in (|40| and (|4T|) . 



Xoe-^(^'^)$(di(x, y)) - K^d2{x, y)) Gt{x, y)dydx, 
K<!>{-d2{x, y)) ~ Xoe-^("'^)$(-rfi(.T, y)) 



Gt{x, y)dydx. 

□ 



4.2.2. Model with BES(3) volatility starting from 1. In this case we formulate a 
condition under which the probabilistic representations hold. 

Proposition 4.4. Assume p <0. If Y is a BES(3) process starting from 1, then 
Xt is a martingale and the probabilistic representations of the arbitrage prices of 
vanilla call and put options hold. 
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Proof. ((2TI) implies that /p Y^dWu is a martingale. Using ([3]), representation 
and 1I35I) we obtain 



Formula (|25|) and p < imply that the local martingale under the last expectation 
is bounded by Xge"^^'^"^"'"^^/^, so it is a martingale and this expectation is equal to 
one. Therefore EXt = Xq and X is a martingale. The assertion of the theorem 
now follows from Theorem 14.11 □ 



5. Appendix 

In this appendix we provide some interesting identities for the Laplace transform 
of [Bt, /g B^du) for a Brownian motion B starting from any point a; S M. The first 
part of Theorem 15.11 is well-known (see Mansuy and Yor [l3l p. 18]). We give here 
a short proof for completeness. The second part is a conditional version of the 
Donati-Martin and Yor formula [3]. Finally, Corollarv 15.21 establishes the density 
function of the vector {Bt, B'^du) for a standard Brownian motion B. 



Theorem 5.1. Let Bt be a Brownian motion starting from x G 
hility space (fi, JF, P) . Then for any c > 0, 6 > and y e M; 



on some proba- 



(42) E^cxp( -cBf-^ f Bldu 



t o / 

^ Jo 

1 



y^cosh(fet) + ^ sinh(te) 



cxp X' 



(43) E 



cxp I - ^ f Bldu 



(6/2 + c)e 



bt 



2 cosh{bt) + f sinh{bt) 



I bt 

sinh(6<) V 2t 



Bt=y 



y^{btcotli{bt) - 1) + 2btx{x + y) 



cosh(&i) — 1 
sinh(6t) 



Proof. Observe that Bt — Vt+x, where Vt is a standard Brownian motion under . 
Let P* be a new probability measure defined by 



(44) 



1^^ = cxp (^bj^ {Vu + x)dVu ^f'^^ 



Then, by Girsanov's theorem, V^ — Vt ~ b JgiVs + x)ds is a standard Brownian 
motion under P*. Moreover, Vt = x[e^^ ^ 1] + e^^ j^e^^'^dV* (see, e.g., Karatzas 
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and Shreve (5). Consequently, changing measure we deduce 



(45) E(^exp cB^ - ^ Bldu 



Ep. exp - c{Vt + xY 



= Ep- ( exp ( - c(yt + xf -b 



(K + x)dV, 
{Vt + xY^t-x^ 



Ep. exp 



exp 



Ep. exp - (&/2 + c)e^^\x + / er'^'dV*) 



cosh(6t) + ^ sinh(6t) 



: exp x 



(6/2 + c)e 



2 cosh(6t) + f sinh(6t) 



as /g e '"'(iV7 is a Gaussian random variable (with respect to P*) with mean and 
variance ^[1 — e"^''*], and for a Gaussian random variable U with mean fi and 



variance a , for any a > 0, we have 



Ep.e" 



: exp 



V2acr2 ^ 1 

So l(42|) is proved. After some algebra we obtain 



2acr2 + 1 / ' 



(46) 



E exp - cB, 



52 /-t 



1 



: exp 



cx + X 



2 6 



(c2 - 674) 



(47) 



y^cosh(6t) + ^ sinh(5t) 

We now prove the second identity. Let us define, for < > 

9 



cot\i{bt) + 2c/b 



Hb,tAy) ■■= E 



exp 



Bf.du 



Bt = y 



It is clear that for i = we have Hbfl.x{y) ~ 1 for any y e M and < Hb,t,x < 1- 
The function Hb,t,x{-) must satisfy, for all c > 0, 



(48) 



E( cxp( -cBt-- / Bidu 



E 



Hb,t,x{Bt) 



It is easy to see that there can only be one such function Hb,t,x- Let us try to find 
Hb t X in the form 



(49) 



Hbj,x{y) = 



F{b,t,x)exp ( - — 



{K{b, t, x) -l){y- xY + 2L{b, t, x){y - x) + R{b, t, x) 



for some functions K, L, R, F of three variables b, t, x. To ease calculations we omit 
the arguments of the functions K, L, R, F. Recall that Bt ~ Vt + x, where Vt is a 
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standard Brownian motion under P, so 

1 



~cB1 



Hb,t,x{Bt) 



F- 



— oo V 27rf 
1 



F cxp — c(x + y) 



1 



Ky' + 2Ly + R] \dy 



cxp 



1 

2i 



2t 

[L + 2ctxf 
K + 2ct 



R 



We now match the last expression with (|46| and guess that 

A' = K{b,t,x) = K{b,t) = btcoth{bt), 
R = 2Lx, 

= - 2LKx + b^xH^. 

The last equation has two solutions Li = xbt{cosh{bt) + l)/sinh(6i) and L2 = 
xbt{cosh{bt) — 1)/ sinh(6t) but the first one tends to infinity when t tends to zero 
(causing to be in the limit). But from the Lebesgue theorem we clearly see 
that 

liuiEHb^tABt) = limE^^cxp - ^ ^ Bldu^^ = 1, 

so what is left is the second solution that fits perfectly our computation. Hence, 
with 

cosh{bt) — 1 



L = L{b,t,x) = xbt- . 

smh[bt) 

R = R{b, t, x) = 2xB{b, t, x) = 2x^bt 



2 cosh(fei) - 1 



sinh(6t) ' 

and F, K defined above, the function Hb,t,x satisfies the condition (|48|) by construc- 
tion. To conclude the proof we insert F, K, L, R just computed in (|49l) and obtain 
the assertion of the theorem. □ 

Corollary 5.2. Fixt > 0. Let g be the density function of the vector {Bt, /J B^du), 
where Bt is a standard Brownian motion. Then for z G M and c > 0, 



(50) 
where 

(51) 



e ''yg(z,y)dy = H*{t,c,z), 



H*{t,c, z) := 



1 



2c 



2tt y sinh(iV2c) 



cxp 



2ccoth(t\/2c; 



Proof. This is an immediate consequence of formula (|43| with x = 0, b 
the fact that 



2c and 



(52) 



/ e-''yg{z,y)dy ^ H^^f^{z)gBt{z), 
Jo 



where gBti^) = ^^^^ denotes the density function of the random variable 



□ 



Bt, and Hbd,x is defined by (|47l) . 
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